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A systematic error in the extraction of sin 2 6w from nuclear deep inelastic scat- 
tering of neutrinos and antineutrinos arises from higher-twist effects arising from 
nuclear shadowing. We explain that these effects cause a correction to the results 
of the recently reported significant deviation from the Standard Model that is po- 
tentially as large as the deviation claimed, and of a sign that cannot be determined 
without an extremely careful study of the data set used to model the input parton 
distribution functions. 

In a recent and stimulating paper [1], the NuTeV collaboration reported a determination of 
sin 2 9w, based on a comparison of charged and neutral current neutrino interactions with a nuclear 
target (Fe), which differs from the Standard Model prediction by three standard deviations. In view 
of the importance of such a result it is vital that the sources of systematic error be clearly identified 
and examined. Here we explain that there is a nuclear correction, arising from the higher-twist 
effects of nuclear shadowing [2,3], for which no allowance has been made in the NuTeV analysis. 
This correction may well be of the same size as the reported deviation. 

The measurement under consideration involves the separate measurements of the ratios of 
neutral current (NC) to charged current (CC) cross sections on Fe for v and v. The best values 
of sin 2 9w and po are extracted from the precisely determined ratios. But the nuclear effects must 
be removed. Because a substantial fraction of the NuTeV data in the region x below 0.1 is at 
relatively low Q 2 (even though the average Q 2 is 16 GeV 2 ), one expects a significant shadowing 
contribution from vector meson dominance (VMD) [4,5], which is of higher twist. As explained by 
Boros et al. [2,3], the effect of the VMD contribution to nuclear shadowing in neutrino interactions 
is substantial, and leads to a reduction of the CC v cross section by about 50% compared with the 
reduction found for photons. (Briefly, the VMD contribution to shadowing is dominated by the p 
meson and f 2 + = 2f 2 , whereas the CC to photon cross sections are in the ratio 18/5.) Together 
with a full NLO analysis of the data, this was important in reconciling the NuTeV and NMC data 
without any need for substantial charge symmetry violation of the parton distributions [6] . A recent 
re-examination of the role of vector meson dominance in nuclear shadowing at low Q 2 finds that 
models (such as that used here) which incorporate both vector meson and partonic mechanisms 
are consistent with both the magnitude and the Q 2 slope of the shadowing data [7]. 

For present purposes we need also to consider this higher-twist effect of shadowing of vector 
mesons for V interactions. These involve predominantly anti-quarks and the shadowing effect is 
relatively larger by a factor of three or so. However, the VMD contribution to shadowing for 
neutral current interactions is 1/2 of that for charged current interactions because Z conversion to 
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a p° occurs with a factor of (1/2) of that for W + — > p + . 

Let us examine how these differences in shadowing effects influence the extraction of sin 2 9 W . 
Suppose the nuclear cross section for NC interactions of neutrinos is larger than that for CC in- 
teractions by a factor of 1 + |e and that the one for anti-neutrinos is larger by a factor of 1 + |e, 
with e expected to be substantially larger than e. Then the nuclear ratios of neutral current (NC) to 
charged current (CC) cross sections are 

where r = a{vN -> £+X)/a(z/iV -> M), r A = r(l + |e)/(l + §e), y£ = 1/2 - sin 2 6 W + 
5/9 sin 4 9 W and ^ = 5/9 sin 4 9 W . Equations (1) and (2) tell us that the nuclear- shadowing cor- 
rections for R A and R^ A are not the same, and that the extraction of sin 2 Q w requires the separate 
knowledge of e and e. 

A detailed analysis of the NuTeV data requires that one model the ratios R A and R^ A at a re- 
quired accuracy of a fraction of a percent. This, in turn, requires an even more accurate knowledge 
of both the quark and antiquark parton distribution functions (pdfs). In general, the pdfs are de- 
rived from a global analysis data from electron/muon, CC neutrino and NC neutrino deep inelastic 
scattering on protons, deuterons and nuclei. The range of Q 2 , particularly at low x (x < 0.1) can be 
quite low. Higher twist shadowing corrections are almost universally ignored in global determina- 
tions of the pdfs. This is certainly the case for the pdfs used by NuTeV. Since the VMD shadowing 
corrections are different for electrons, CC neutrino scattering and NC neutrino scattering, the pdfs 
resulting from such a global analysis are at best an approximation to the true ones, with unknown 
systematic errors. Worse, one cannot simply add a shadowing correction to a simulation based on 
such global pdfs as even the sign of the correction will depend on the particular data sets included 
in the analysis. Of course, the systematic errors encountered will not be serious for most purposes. 
However, in this case, where the signal of a deviation from the Standard Model is at the percent 
level, one must control this potential source of error extremely carefully. 

We cannot undo the global analysis of pdfs by the NuTeV collaboration. However, we can 
make an estimate of the order of magnitude of the shadowing corrections using Eqs. (1) and (2). 
The quantity e is given as a function of x for Q 2 = 5 GeV 2 as the dashed curve of Fig 3b. in 
Ref. [2]. One needs to take the deviation between dashed curve and unity. Thus, e.g. e(x, Q 2 = 
5GeV 2 ) 0.041 at x = 10~ 2 . For other values of Q 2 , one may use: e(x, Q 2 ) ps e(x, Q 2 = 

5GeV 2 ) ( m j;t+Q C 2 V2 ) • Furthermore e/e ^ F^{v)/F^{v) ps 2. Using a typical value of Q 2 ps 

lOGeV 2 and range of x ps 0.05 of the NuTeV experiment leads to the values e = 0.006 and 
e = 0.012 [8]. Thus the nuclear corrected ratios R A and R A would be smaller than those reported 
in Ref. [1], by 0.003 and 0.006 respectively. That these numbers represent important corrections 
can be seen immediately by examining the sources of errors reported in Ref. [1]. The total error 
for R\ is reported as 0.0013 and that for R A as 0.00272. The effects of shadowing are larger than 
these quoted errors by a factor of two or three! It is clear that any analysis of nuclear data aimed at 
determining sin 2 9 W must account for nuclear shadowing. But this has not been done in Ref. [1]. 

It is necessary to carry out a more refined analysis of the data which properly incorporates 
the high twist components of nuclear shadowing, starting with the pdfs themselves. Such an 
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analysis should also take into account the experimental acceptances as a function of x, y and Q 2 . 
Alternatively, one could drastically reduce the VMD contribution by restricting the data set to 
events with Q 2 > 5 GeV 2 - although we understand that this may present difficulties for NC 
events. 

We find that the size of the shadowing effects is substantial and should be incorporated in the 
experimental analysis. It is true that the simplest estimate of this effect (ignoring the effects on the 
pdfs themselves, which was discussed above) is opposite to that required to explain the deviation 
from the Standard Model. Thus it may be that the deviation from the Standard Model could be 
even larger than reported in Ref. [1]. 

Finally, we note that several other effects that tend to reduce the discrepancy have been re- 
ported. The influence of charge symmetry breaking, arising from the mass difference between up 
and down quarks [9], accounts for about a third to a half of the deviation between NuTeV's value 
of sin 2 ^ and that of the Standard Model in a model-independent manner [10] . Furthermore, 
it has been known for more than 20 years that parton distributions of nucleons bound in nuclear 
matter differ from those of free nucleons. Such effects [11] still present a considerable challenge 
to our understanding of nonperturbative QCD and it is not inconceivable that they could eventually 
account for the entire deviation of sin 2 9w- 

It seems clear that the extraction of the the value of sin 2 9 W from neutrino-nuclear interactions 
involves handling several different types of corrections of different signs, including some that 
are difficult evaluate with precision. The situation here may well be similar to many in strong 
interaction physics, in which a "cocktail" of effects is required [12]. Considering that possible 
explanations in terms of new physics are not compelling [13], considerable efforts must be applied 
before concluding that the NuTeV result really demonstrates a deficiency of the Standard Model. 
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